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Abstract 

Within the "pseudo-Dirac" scenario for massive neutrinos the existence of sterile 
neutrinos which are almost degenerate in mass with the active ones is hypothesized. 
The presence of these sterile neutrinos can affect the flavor composition of cosmic 
neutrinos arriving at Earth after traveling large distances from astrophysical objects. 
We examine the prospects of neutrino telescopes such as IceCube to probe the very 
tiny mass squared differences 10~^^ eV^ < Am^ < 10"^^ eV^, by analyzing the ratio 
of )U-track events to shower-like events. Considering various sources of uncertainties 
which enter this analysis, we examine the capability of neutrino telescopes to verify the 
validity of the pseudo-Dirac neutrino scenario and especially to discriminate it from 
the conventional scenario with no sterile neutrino. We also discuss the robustness of 
our results with respect to the uncertainties in the initial flavor ratio of neutrinos at 
the source. 

PACS numbers: 14.60.Pq; 13.15.+g; 95.85.Ry 
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1 Introduction 



Analyses of the data from reactor [T], accelerator p], atmospheric [3] and solar [1] neutrino 
experiments conclusively demonstrate the oscillation of neutrino flavors. The results of these 
experiments can be interpreted by two independent mass squared differences (between three 
active neutrinos). From the direct measurement of the invisible part of the decay width of Z 
boson (z.e., Tz^u^uc,), the number of active neutrinos lighter than Mz/2 found to be N^, = 
2.92 ±0.06 [HI [6] and from the fit of the LEP data to Standard Model prediction it found to be 
Nu = 2.994±0.012 [3 [6]. Thus, if an extra light neutrino exists, it should be a sterile neutrino 
(singlet under the gauge symmetries of the Standard Model.) Historically the strongest hint 
for the existence of sterile neutrino came from the short baseline LSND experiment [8]. Data 
of the LSND experiment, with neutrino energy ~ 30 MeV and baseline ~ 30 m, suggested 
a 3+1-scheme (active+sterile) with the new mass squared difference Am^ ~ C^(l) gV^. The 
LSND result has not been verified by the MiniBooNE experiment P| and considerable efforts 
have gone into reconciling the null result of MiniBooNE with the LSND data [10] . All the 
other data of the neutrino experiments can be interpreted by assuming only three active 
massive neutrinos without any need to introduce sterile neutrinos in the data analyses. 
However, sterile neutrinos can still be present in the yet not probed regions of the parameter 
space (Am^, 6). These regions correspond to sterile neutrinos almost degenerate in mass with 
the active ones, with very tiny mass differences Am^ -C Am^g^. The scenario of degenerate 
sterile neutrinos, the so-called "Pseudo-Dirac" % scenario, has been proposed long time ago 
in [11] and has been studied in the literature extensively |12j . 

The prospect for the existence of light sterile neutrinos with masses nearly degener- 
ate with the masses of active neutrinos is motivated in many theoretical extensions of the 
Standard Model [I3j. From the observational point of view, probing very small Am^ be- 
tween sterile and active neutrinos needs very long baselines. Neutrinos coming from the 
Sun (which is the farthest observed source of neutrinos with continuous emission) set the 
bound Am^ < 10^^^ eV^ on the active-sterile mass splitting. Bounds from other performed 
or forthcoming experiments will be discussed in Sect. 12.11 In this paper we concentrate on 
the effects of almost degenerate sterile neutrinos on the expected flux of cosmic neutrinos 
coming from astrophysical sources. 

The new generation of km^ scale neutrino telescopes give a unique opportunity to probe 
the very tiny Am^ in pseudo-Dirac scenario. The cosmic neutrinos from sources such as 
GRBs [E], AGN [15] and type Ib/c supernovae [16] travel large distances over ~ 100 Mpc 

^The reason for this nomenclature will be described in Sect. [5] 
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before arriving at neutrino telescopes in Earth. With such extremely long baseline, tiny 
mass squared differences as small as Am^ ~ 10^^^ eV^ {E,^/1QQ GeV) can be probed. The 
idea of using neutrino telescopes to discover the sterile neutrinos present in pseudo-Dirac 
scenario was proposed in [TTl [HI [191 (SO]- In order to probe small values of Am^, it has 
been suggested to look at distortions in the spectrum of from supernovae remnants in 
the average distance of ~ 1 — 8 kpc in [T7] and the spectrum of from Galactic center in 
[18] . The authors of [HI [20] evaluate the effect of the pseudo-Dirac neutrinos on the flavor 
composition of the cosmic neutrinos; i.e., the deviation of the F^^ : F^^ : F^,^ (where F^^ 
is the flux of Ua + i^a at Earth) from the expected value 1 : 1 : 1 in the absence of sterile 
neutrinos. 

In the measurement of flavor composition of neutrinos in neutrino telescopes, the un- 
certainties in the relevant parameters and experimental limitations should be taken into 
account. For example, there are uncertainties in the mixing parameters of neutrinos and 
also in the spectrum of the arriving neutrinos. Also, the current constructed or proposed 
neutrino telescopes, AMANDA/IceCube [21], NEMO [22], NESTOR [23], ANTARES [24j 
and KM3NET [25] cannot identify all three flavors of the active neutrinos. Ref. [26] consid- 
ers these uncertainties and experimental limitations in the analysis of the cosmic neutrinos 
in order to extract mixing parameters and flavor composition of neutrinos at the source. In 
this paper, by considering the aforementioned uncertainties and experimental limitations, 
we investigate the potential of neutrino telescopes in discovering the pseudo-Dirac nature of 
neutrinos. The initial flavor ratio of the neutrinos at the source can also be a source of un- 
certainty in the calculation of event rates in neutrino telescopes. We discuss the robustness 
of our result to this kind of uncertainty. 

The paper is organized as follows. In sect. [2], the pseudo-Dirac scenario for massive neutri- 
nos is reviewed. In sect. 12.11 the current bounds on Am^ from various neutrino experiments 
are summarized; and in sect. 12.21 the effects of sterile neutrinos on the flavor composition of 
cosmic neutrinos are discussed. Sect. [3] is devoted to the production mechanism of neutrinos 
at the source and their detection processes in the neutrino telescopes. Various sources of 
uncertainties that enter the calculation of event rates in neutrino telescopes are enumerated. 
SectJH summarizes the results of the present analysis on the capability of neutrino telescopes 
to discriminate between pseudo-Dirac and conventional scenarios. A summary of the results 
and the conclusions are given in sect. [5l 
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2 Pseudo-Dirac Scenario 



A simple and economic way to generate mass for neutrinos in the SM is to add right-handed 
(sterile) neutrinos to the matter content of SM. In the presence of A^^ right-handed (sterile) 
fields u^R {k = 1, . . . , Ng), we define the following column matrix \E' of = 3+A^s left-handed 
fields 

^ = (z^eL, Z^mL, T^tL, (z^Ir)^, • • • , (z^TV.i?)^)^ , (1) 

where u'" = Cv^ and C is the charge conjugation operator. For Majorana neutrinos which 
we consider here {vaR)'" = i^aL- Here we consider models with at most three sterile neutrinos 
[Ng < 3). In the basis the generic mass term for neutrinos is 

£^ = -^^M^ + H.c., (2) 

The (3 + A^s) X (3 + A^s) mass matrix M is of the following form (after electroweak symmetry 
breaking) 

M = ( < ] . (3) 

where is the A^^ x 3 Dirac mass matrix and and mpi are the 3x3 left-handed and 
Ng X Ng right-handed Majorana mass matrices, respectively. The non-vanishing elements 
of mi and uir violate lepton numbers while by assigning the lepton number -|-1 to sterile 
neutrinos, mo conserves this symmetry. The left-handed mass matrix is not invariant 
under the SM gauge group SU{2)l and should be zero unless other new particles (such as 
a new Higgs triplet) are present. The elements of can take a wide range of values, it 
can be as large as the GUT scale ~ 10^^ GeV which are preferred in see-saw mechanisms, 
or it can vanish like tjil as a result of new gauge symmetries such as SU{2)ji [27]. The 
case rriL = uir = and A^^ = 3 results in pure Dirac neutrinos. In this case the six Weyl 
neutrinos decompose into three pairs of neutrinos with degenerate masses. The active-sterile 
mixing angle in each pair is maximal 6 = 7r/4, but the active neutrinos do not oscillate to 
their sterile partners because Ami n = ^1 —^"i = 0, where mo. and mn are the masses of 
sterile and active neutrinos in the j-th pair, respectively. Here we are interested in the case 
^Li^R m£,. The non-zero but very small values of the elements of and lift the 
degeneracy in mass at each pair. In this "pseudo-Dirac" scenario, active-sterile mixing angle 
in each pair is ^ ^ 7r/4 and active-sterile oscillation can in principle occur due to very small 
but non-zero Am^^. To illustrate this point, let us consider the one generation example. In 
this case, the mass matrices m^, m^ and in Eq. ([3]) are numbers (we assume that all 
the masses are real.) In the pseudo-Dirac limit, we obtain tan(26') = \2mo/ {mpt — m^)! ^ 1 
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and Am^^ ~ 2mj:){mL + m/j) <^ m^- Notice that in the pseudo-Dirac scenario, neutrinos 
oscillate even in one generation, in contrast to pure Dirac scenario where oscillation occurs 
only between generations. 

In general the N x N symmetric mass matrix M (where N = S + Ng) can be diagonalized 
by VjMVy = Mdiag, where K is a x unitary matrix. We choose the elements of 



diag(ma^ , rua^ , , , 



m 



The mixing matrix V appearing in 



such that Mdiag 

the weak charged-current J^r = 2'$iV^'-f'^laL is a 3 x rectangular matrix with the elements 
Vak = l^/3=e /X r(^z^)a/3(^i')/3fc; whcrc Vi is the 3x3 diagonalizing unitary matrix of charged 
leptons mass matrix. In the case A^^ = 3, the 6x6 matrix V can be parameterized by 12 
mixing angles and 12 CP- violating phases (7 Dirac phases+5 Majorana phases.) It has been 
shown in [28] that in the pseudo-Dirac limit m^, m^j <^ rriu and at first order of perturbation 
in the small parameters niL/mo and m^j/m^), the mixing matrix V has only three mixing 
angles (responsible for oscillation between the pairs) and three CP- violating phases (1 Dirac 
phases-|-2 Majorana phases). This fact can be seen from the explicit form of the matrix Vy 
which diagonalizes the mass matrix M (in the pseudo-Dirac limit) [28] : 



PMNS 









1 



'3x3 



3xNs 



u 



R 



'V2 I (X 



3xNs 



-il 



(4) 



NsXNs 



where Upmns is the 3x3 conventional neutrino mixing matrix of left-handed neutrinos, Upt 
is the Ns X Ns unitary matrix which diagonalizes the right-handed Majorana mass matrix, 
Inxn is the n X n identity matrix and the matrices Xsxat^ {Ns < 3) are: 



^3x1 



/ 1 \ 





3x2 



/l o\ 

1 

Vo o; 



X3x3 — -^3x3- 



(5) 



The flavor conversion probability between the active neutrinos Pa/3 = Pua^u^iL, E^) is 



a(3 



Kexp 



2E„ 



al3 



(6) 



Using the explicit form of the matrix Vy in Eq. (jl]), the probability P^/? becomes 



P 



(7) 
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where and mj are the mass eigenvalues in the j-th pair of active and sterile neutrinos; 
Uaj and U/3j are the elements of the 3x3 mxixng matrix Upmns- Notice that this relation 
reduces to the standard flavor conversion probability formula in the limit of pure Dirac 
neutrinos = mJ (j = 1,2,3). By setting = mJ for the active neutrino generations 
which do not have sterile partners, Eq. ([7]) also applies to cases with Ng < 3. 

Using Eq. ([7]) in analyzing the data of oscillation experiments gives information on 
{m^)'^ — {'m~Y in each pair. In subsect. 12.11 we review the current bounds on active-sterile 
mass square differences and the prospect of future experiments to improve these bounds. 
In subsect. 12.21 we discuss the implications of Eq. ([71) on the flavor composition of cosmic 
neutrinos. 

2.1 Current Bounds on Am^ and Sensitivity of Future Experi- 
ments 

An oscillation experiment with baseline L and neutrino energy Ey can probe mass square 
difference Am^ ~ E^/lAirL). If Am^ <^ Ey/{ATTL), the baseline is too short for flavor 
oscillation to take place; on the other hand, if Am? ^ Ey/ (AirL) so many oscillations take 
place during the propagation and the oscillatory term should be averaged out. In both of 
these cases it is not possible to derive the value of Am^ in oscillation experiments. 

Solar neutrino experiments with the baseline 1 AU ^ 1.5 x 10^^ m and neutrino energy 
Ey ~ 0.1 — 10 MeV, can probe mass squared differences Am^ ~ 10^^° — 10~^^ eV^. These very 
small values of Am^ has been favored by the so-called "Vacuum Oscillation Solution" of the 
solar neutrino problem, but as is well-known this solution has been ruled out by KamLand 
[1]. However, the sterile-active oscillation with mass square differences Am^ < 10^^^ eV^ can 
still be present as a subdominant effect in solar data. The recent work [29] updates the solar 
data and obtains Am^ < 1.8 x 10~^^ eV^ (at 3a level) for the sterile-active mass splitting. 
This bound is the most stringent bound on Am^. The flavor composition of the neutrinos 
from core-collapse supernovae (SNe) also can change due to an active-sterile oscillation from 
the SN to Earth. The mean energy of the neutrinos from a SN explosion is E^ ~ 30 MeV. 
Thus, a SN explosion at a distance of ~ 10 kpc can probe Am^ ~ 10~^^ eV^. The constraint 
from the data of the SN1987A data is not restrictive because of the low statistics and high 
uncertainties in the mechanism of SNe explosion [30]. Construction of future Mton water- 
Cerenkov detectors can dramatically improve the current bound or find a hint for sterile 
neutrinos hypothesizes in pseudo-Dirac scenario [31] . 
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Population of the sterile neutrinos in the early universe and their effects on the Big 
Bang Nucleosynthesis (BBN) can change the abundance of light elements. The standard 
BBN, given the number of the relativistic particles N^, and the baryon asymmetry r] = 
ns/n^, predicts the abundance of the light nuclei in the universe. Assuming that the sterile 
neutrinos are produced only in the active/sterile oscillation and the initial abundance of 
sterile neutrinos at temperatures T 3> MeV is zero, the tightest limit comes from the '^He 
abundance: Am^ < IQ-^ eV^ [3ll[32]. 

Two main non-oscillation neutrino experiments which probe neutrino masses kinemati- 
cally are tritium beta decay and neutrinoless double beta decay (Oz//5/5) experiments. Among 
them, the OujSjS decay is sensitive to the Majorana or Dirac nature of neutrinos. The rate of 
0uj3f3 decay is proportional to the effective mass of the electron neutrino which is defined as 



As mentioned after Eq. ([3]), in the limit rriL = fn^ = (pure Dirac neutrino) each Dirac neu- 
trino is the superposition of two Majorana neutrinos with degenerate masses and opposite 
CP eigenvalues. It is easy to see that the Majorana neutrinos in each pair interfere destruc- 



pseudo-Dirac scenario with non-zero Majorana masses and mL^mn -C m^i, the cancelation 
is not exact and {rriee) 7^ but it is very small |_33j. Thus, it seems that the observation of a 
positive signal in the next generation 01/(3(3 experiments, with sensitivities (mee) ~ 10 meV, 
will rule out the small values of and rriR in the mass matrix of neutrinos and therefore 
pseudo-Dirac scenario. However, two points should be considered. The first one is that the 
value of (mee) can still be significant if only one or two families of neutrinos have sterile part- 
ners. Contribution of each family to the value of (mee) depends on the corresponding mixing 
matrix element Uaj- This means that, because of the small value of f/es 0.041), presence 
or absence of a sterile neutrino with a mass almost degenerate with u^l do not change the 
value of (rriee) substantially; but the case with two sterile neutrinos with a masses degenerate 
with I'll and z/2l leads to an effective mass (mee) much smaller than its value in the absence 
of sterile neutrinos. The second point is that the dominant contribution to 0^(3 (3 decay can 
come from new particles or physics beyond the SM, such as a. V + A interaction [3lj. It is 
shown in [35] that a non-zero Oi>(3f3 decay rate generates small m^ through radiative correc- 
tions, which results in pseudo-Dirac scenario for neutrino masses. Considering these points, 
it is not easy to draw a conclusion on pseudo-Dirac scenario from the results of the OujSP 
experiments. 
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tively in Eq. ([8]) {rrij 



+ 



rrij ) which results in (mee) = for pure Dirac neutrinos. In the 
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2.2 Cosmic Neutrinos 



Neutrinos arriving at neutrino telescopes from astropliysical sources travel distances of the 
order L ~ 100 Mpc. The flavor conversion probabilities over these large distances can be 
obtained by averaging out the oscillatory terms in Eq. ([7j). Two different scales of Am^ are 
involved in Eq. ([7!), one is the atmospheric, Am^j^ ~ 10~^ eV^, and solar, Am^^^ ~ 10~^ eV^, 
and the other is the very small Am? between the mass eigenstates in each pair of active and 
sterile neutrinos. As it is shown in [36] the oscillatory terms depending on Am^^^ and Am^^; 
should be completely averaged out over these large distances. The mass squared difference 
Am? that can be probed by neutrinos with energy which propagate through distance L 
is 

Thus, even for L as large as 10 Mpc oscillatory terms given by Au? ~ 10^^'' eV^ do not 
average out. After averaging out Am\^^ and Am^^;, Eq. ([7]) becomes [19] 



Y.\U^,?\U,,?cos'[-^y (10) 



a/3 — 



where Am| = {rn^Y ~ ("^/^^ mass squared difference in the j-th pair. Thus, if the 

initial flavor composition of neutrinos in the source is We '■ Wf^ : Wr, the flavor composition of 
the neutrino beam arriving at Earth will be F^,^ : F^^ : F^,^, where 



/3 i=i 

The average of the cosine factor for Am,'jL/4E^ ^ 1 is 1/2. Thus, if for all three pairs 
(j = 1, 2, 3) the condition Am'^L/AE^, ^ 1 is satisfied, all three F^,^ in Eq. ffTTl) are multiplied 
by 1/2 such that the flavor ratios F^^ : F^^ : F^^ do not change with respect to the flavor 
ratios in the pure Dirac case Am?- =0. In this situation the only difference between the 
pure Dirac and pseudo-Dirac scenarios is that the number of neutrinos arriving at Earth 
is reduced by half in the pseudo-Dirac case. Thus, it is very hard to verify pseudo-Dirac 
scenario for distances or mass squared differences where Am'^L/AEu ^ 1 for j = 1, 2, 3. That 
is because the estimation of the overall normalization of neutrino flux needs knowledge about 
the details of the neutrino production mechanism in the source which is not well understood. 



8 



However, the pseudo-Dirac scenario of neutrinos can be tested in neutrino telescopes for 
the cases that only one or two sterile neutrinos exist (A^^ < 3) or the distance of the source or 
the mass square differences are such that all the three oscillatory terms given by ArrijL/AEi, 
do not average out. In these cases the average of one or two of the cosine factors in Eq. (fTT!) 
are 1/2 and the other cosine factors can be replaced by 1 (we do not consider the special 
situations where Am'jL/AE,^ ~ 1). Thus, measuring the deviations of the flavor ratios in 
Eq. (fTT]) from their standard values (in the absence of nearly degenerate sterile neutrinos) in 
neutrino telescopes can shed light on these cases. In the next section we discuss the details 
of the detection processes in the neutrino telescopes and the feasibility of identifying differ- 
ent neutrino flavors in these experiments. By taking into account the realistic measurable 
quantities in neutrino telescopes and the uncertainties in these measurements, we discuss to 
what extent it is possible to measure the flavor ratio of cosmic neutrinos and their deviations 
from the standard values. 



3 Detection and Production Processes 

In this section we briefly discuss neutrino flavor identification in the km^ scale neutrino 
telescopes such as IceCube or its counterparts in the Mediterranean sea; KM3NET, NEMO, 
NESTOR and ANTARES. A detailed description of the flavor tagging efficiencies in a typical 
neutrino telescope can be found in [261138] . Here we summarize the main points relevant for 
the present analysis. Particularly, for the first time, we take into account different sources 
of uncertainties in the calculation of event rates and also a more realistic analysis of the 
detectable events, such as the contribution of z/,- (z/^) to the yU-tracks. 

The flux of neutrinos arriving at Earth can come from a single luminous point source 
or as a diffuse flux from sum over different sources at different distances. The advantage 
of the point sources is that the short period of the burst and the direction of incoming 
neutrinos can be used to reduce the background events (especially when the source can be 
identified using a different method such as gamma photons for GRBs). Point sources with 
an intense neutrino flux detectable at km'^-scale neutrino telescopes can take place in the 
close-by galaxies located at a distance of < 10 Mpc and such a source of neutrinos yields 
about a few hundred neutrino events in IceCube [37] . 

Discriminating between different flavors of neutrinos is a great challenge for neutrino 
telescopes. Two types of events are completely distinguishable in the next generation of these 
experiments: /x-track events and shower-like events. Charged Current (CC) and Neutral 
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Current (NC) interactions of different flavors can contribute to each of these events, /i-track 
events, which are the Cerenkov hght radiated by muons propagating through the volume of 
the detector, get contributions from two sources: i) /i (/i) produced in the CC interaction of 
ip^); ii) CC interaction of (z/^) which produce r (r) leptons and the subsequent leptonic 
decay of tau leptons r fJ^i^^i^r (^^ U^-'^i^t) produce /i (/i). Shower-like events have three 
sources: i) NC interactions of all the three flavors of neutrinos; ii) CC interactions of z/g 
and Pe] and iii) CC interaction of z/,- {ut) which produce r (r) leptons and their subsequent 
hadronic decays. The exact formulae for calculating the rate for each of these events can be 
found in Sect. 2 of p6] . 

The threshold energy of the detection of /i-tracks and showers in experiments such as 
IceCube, respectively, is Ejj" ~ 100 GeV and ^*L«,er ~ 1 TeV [21]. On the other hand, 
the mean free path of the neutrinos with energy i?^"* ~ 100 TeV becomes of the order of 
~ 2i?0, the diameter of Earth. The exact values of the Ejj^, Ef^^^^^ and i?^"* depend on the 
details of the experiments, such as the geometry of the photomultipliers in the volume of the 
detector and the direction of the incoming neutrinos; and a dedicated analysis can be done 
for each experiment. Here, in order to avoid considering the absorption of neutrinos in the 
Earth, we restrict the analysis to 100 GeV < E^ < 100 TeV. 

The realistic quantity that can be measured in neutrino telescopes is 



The value of R can be calculated if we know the initial flux and the flavor ratio of neutrinos 
at the source. The main mechanism of neutrino production at the astrophysical sources is 
the interaction of the accelerated proton by the ambient protons and photons. The decay of 
the secondary particles (vr^, K^, D, . . . ) produced in these pp and py interactions generate 



The exact value of the flavor ratio We : : Wr depends on the spectrum of the parent 
particles and the properties of the production medium. A class of models based on the 
Fermi acceleration mechanism for the particle in the source, predict a power-law spectrum 
for neutrinos 



R 



Number of Muon-track events 



(12) 



Number of Shower-like events 



neutrinos and muons. For example, the pion chain decays 7r+ f^^i^fi {G^i^ei^fi)i^fi and 
TT^ fi^iy^ — > (e~z>et'^)^'^ generate neutrinos with the flavor ratio (z/^ + z/^)/(z/e + z^e) — 2. 




(13) 



dE, 
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where a is the spectral index and Mvf^ is a normahzation factor. Acceleration of particles 
through Fermi acceleration mechanism [39| results in a = 2 for neutrino spectrum. However, 
non-linear effects change this value such that a can take any value in the interval (1, 3) [40j. 
It is shown in [?T] that for the case of pion decay chain and assuming a = 2, the initial 
flavor ratio is : : = 1 : 1-85 : (the difference with 1:2:0 comes from the wrong 
polarization states of /i^ in the decay of vr^). Also the authors of [12] show that inclusion 
of other secondary particles (such as has a very little effect on this value. The //^ 
generated in the decays of the secondary particles can substantially lose their energy before 
decay. In this case the neutrinos generated in the decay of muons do not contribute to the 
flux of neutrinos with 100 GeV < -Ejy < 100 TeV and the flavor ratio becomes 0:1:0. 

In the calculation of R for different scenarios of neutrino production at source and prop- 
agation between source and Earth, the uncertainties in the input parameters should be 
considered. Uncertainties of the input parameters in the calculation of the /i-track and 
shower-like event rates induce uncertainties in the value of R. Here we summarize these 
sources of uncertainties: 



Mixing Parameters Flavor content of the neutrino beam arriving at Earth depends on 
the mixing angles [Qvi-, ^^23, ^13) and Dirac CP- violating phase b through the | f^aj P | t^/jj P 
factors in Eq. ( ITTi) . An update on the values of these parameters and uncertainties in 
each of them can be found in [46j, which are also listed in Tabled] 

Spectral Index As it is mentioned after Eq. (fT3l) . the spectral index a can take any value 
in the interval (1, 3). It is shown in [38] that IceCube can measure the spectral index a 
with 10 % precision (assuming EldF^j dEy = 0.25 GeV cm~^ sr~^ yr~^ and after one 
year of data-taking). 

i^a/^a Ratio Let us define A/3 = ATp^/AC^ for [3 = e, fi. It is obvious that = 1 in the pion 
decay chain, but the value of Ae depends on the ratio of /n^ production and can take 
any value in the interval (0, 1). To best of our knowledge, in the energy range we are 
interested here (100 GeV < Ey < 100 TeV), there is not any proposed or established 
method to measure Ag. 

Neutrino-Nucleon Cross Section The rate of the /i-track and shower-like events depends 
on the CC and NC neutrino-nucleon cross sections a^^ and cr^^ . The current uncer- 
tainty in these cross sections is ~ 3 %. However, the uncertainties in cr^ and a^^ 
have a very small effect on R because of the cancelation between numerator and de- 
nominator of Eq. f|T2|) . 
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Parameter 


Best-fit 


Current Allowed Range 


Future Uncertainty 


sin^ 612 


0.304 


0.25 - 0.37 (3(T C.L.) 


6 % (Ref. [13]) 


sin^ ^23 


0.50 


0.36 - 0.67 (3(T C.L.) 


6 % (Ref. m]) 


sin^ 6^13 


0.01 


< 0.056 (3(T C.L.) 


5 % (Ref. [15]) 


5 




[0,27r) 


[0, 27r) 


a 




10 % 


10 % 






[0, 1] 


[0,1] 



Table 1: Relevant parameters in the calculation of R with their uncertainties. The current 
uncertainty column represents the 3(T uncertainty interval for mixing angles [16]. The future 
uncertainty column shows the precisions that can be achieved in the forthcoming neutrino 
oscillation experiments, as described in the corresponding references. For ^13 it is assumed 
that its value is close to the current upper limit: sin^ ^13 = 0.03. 



In Table [T] we have listed the relevant parameters in the calculation of R and their 
uncertainty intervals. The third and forth columns respectively correspond to the current 
and future uncertainties in these parameters. In addition to the uncertainties in the input 
parameters listed above, measurement of R is also done with limited precision and has an 
uncertainty. It is shown in [38] that by assuming that the flux of neutrinos is E^dF^/dEi, = 
0.25 GeV cm~^ sr~^ yr~^, the ratio R can be measured with ~ 7 % precision after a couple 
of years of data-taking. 

Now, by taking into account the uncertainties mentioned in this section, the question is to 
what extent it is possible to measure the deviations of the flavor composition F^^ : F^^ : F^^ 
from its value in the absence of pseudo-Dirac sterile neutrinos. In the next section we will 
discuss the prospect of neutrino telescopes to measure these deviations. 

4 Results and Discussion 

In this section we consider two sources with different initial flavor compositions: the pion and 
stopped-muon sources with flavor compositions 1 : 1.85 : and 0:1:0, respectively. The 
value of R in Eq. (1121) in the absence of almost degenerate sterile neutrinos and assuming the 
best-flt values for the mixing angles, 6 = and Ae = 1, will be denoted by -R,r and R^, for 
pion and stopped-muon sources respectively. We assume the power-law spectrum Eq. f|T3l) 
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for the neutrino production, with the spectral index a = 2. At the end, we investigate 
the robustness of the resuhs with respect to the deviation of initial flavor ratios from the 
assumed values 1 : 1.85 : and 0:1:0. 

Fig. ([1]) shows R versus a for pion source, where R^^ = 2.50. The hatched area between 
the two red curves corresponds to the values that R can take if the input parameters vary in 
the uncertainty intervals. The inputs have been varied in the current uncertainty intervals 
shown in the third column of Table [H The two horizontal and vertical pairs of dashed-lines 
show the 7 % and 10 % precision intervals in the measurement of R and a, respectively. 
The rectangle created from the intersection of these dashed-lines corresponds to the region 
of parameter space {R, a) where can be limited by the measurements in neutrino telescopes. 
Thus, the points inside the rectangle represent the values of R for the case of no sterile 
neutrino and consistent with 7 % (10 %) precision in the measurement of R (a). The green 
dashed-curves in Fig. ([1]) show the values of R in the presence of sterile neutrinos. For 
example, in Fig. ([U-a) the hatched region between the green dashed-curves corresponds to 
the case when the sterile neutrino mass is almost degenerate with the mass of the active 
neutrino Uil such that the oscillatory term depending on Amf in Eq. ([71) can be averaged 
to 1/2 and the flavor conversion probabilities become 

Pal3 = l\Ual\^\UfS,\^ + \Ua2\^\Uf,2\^ + | f/„3 l' | f//33 1 ' • (14) 

As can be seen, in all six parts of the Fig. ([1]) the two hatched areas overlap. The overlap- 
ping of the hatched areas inside the dashed-line rectangles means that the presence of the 
sterile neutrinos cannot be ruled out, even if the measurement of R gives a value inside this 
rectangles. However, if the measurement of R gives a value much different than R^^, the 
existence of sterile neutrinos is favored. For example, the value i? = 2.1 is not possible in 
the no sterile case and this value is favored by the scenarios depicted in Fig. ([T]-b,c,f). 

Fig. ([2]) is the same as Fig. with the exception that the input parameters have been 
varied in the future uncertainty intervals shown in the forth column of Table [H We assumed 
(6%, 6%, 5%) uncertainties for (sin^ 6*12, sin^ 6^23, sin^ 6*13), which can be achieved in the forth- 
coming neutrino oscillation experiments [lUllllHS]. Also, we assumed a large value for the 
13-mixing angle sin^ ^^13 = 0.03 which is near its present upper bound. The uncertainties 
of 6, Ae, a and R are the same as in Fig. ([T]). As can be seen from Fig. ([2]), reducing the 
uncertainties of mixing angles results in a separation between the two hatched areas where 
there was an overlap in Fig. ([1]). For example, in Fig. ([2]-f), the regions corresponding to 
no sterile neutrino case and the case with averaged second and third pairs are completely 
separated. This separation means that if the measurement of R gives a value inside the 
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Figure 1: The dependence of R, ratio of ^-tracks to shower-like events, on a for the pion source 
with the initial flavor ratio 1 : 1.85 : and power-law spectrum with the spectral index a = 2. 
Assuming the best-fit values for the mixing angles, 6 = and Ag = 1, the value of R is R-,^ = 2.50. 
In each figure the red curves represent the case with no sterile neutrino and the green dashed-curves 
correspond to the cases with average on pairs mentioned in the legends. The hatched areas show 
the values that R can take when the input parameters vary in the current uncertainty intervals in 
Table [TJ The two vertical and horizontal dashed-lines show the 10 % and 7 % precisions in the 
measurements of a and R, respectively. 
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Figure 2: The same as Fig. ([T]) except that the input parameters have been varied in the future 
uncertainty intervals of Tabled) Particularly, we assumed 5 % uncertainty interval for the 13-mixing 
angle with the central value near the present upper bound: sin^ ^13 = 0.03. 
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rectangle in Fig. ([2]-f), the existence of sterile neutrinos (for the case of average on second 
and third pairs) will be ruled out. On the other hand, a value of R outside the rectangle 
can be interpreted as a signal for the existence of sterile neutrinos. However, recognizing 
which case is consistent with the measurement of R should be done with care. To illuminate 
this point, let us assume that IceCube has measured R = 3.0. According to the diagrams in 
Fig. ([2]) this value of R is consistent with three cases: i) average on first pair; ii) average on 
first and second pairs; and iii) average on first and third pairs. Notice that with the current 
uncertainties of the mixing angles (depicted in Fig. ([1])) it is not possible to conclude that 
R = 3.0 is a signal of sterile neutrinos. 

Figs. ([3]) and (jlj) show the dependence of i? on a for the stopped-muon source; where 
R^ = 3.07. The input parameters have been varied in the current uncertainty intervals 
in Fig. ([3]). Comparing Fig. ([3]) with Fig. ([T]), it is obvious that the hatched areas are 
wider for the stopped-muon source, which means that the recognition of sterile neutrinos is 
harder for this kind of sources. In contrast to the pion source, reducing the uncertainty of 
input parameters to the future uncertainty intervals do not lead to a separation between the 
hatched areas in Fig. (jll). However, some cases can marginally be discriminated such as the 
case shown in Fig. (jU-d). 

In drawing Figs. ( |T|2|) and Figs. ( 13|4|) we have assumed initial flavor ratios 1 : 1.85 : and 
0:1:0, respectively. However, it should be noticed that the initial flavor ratio can deviates 
from these values due to the interplay of different mechanisms in neutrino production or 
effects that have not been considered in the calculation of these values. For example, in the 
pion source, a part of the produced muons in the decays of tt^ (not all of them) can lose their 
energy before decay such that the initial flavor ratio of neutrinos takes a value between the 
two extreme cases of pion and stopped-muon sources. The question that arises here is that 
to what extent the results of this section (for example the separation between hatched areas 
in Fig. (EJ-f)) are robust against the deviations of the initial flavor ratio We '. : Wr from 
the assumed values To answer this question we consider particularly the case of Fig. ([2]-f) 
where the hatched areas are completely separated. Also we use the following parametrization 
of the initial flavor ratio 

We : Wfj, : Wr = n : 1 : 0. (15) 

The Wr = in the above parametrization comes from the fact that the number of prompt 
Ur (from the decays Dg —>■ rur, . . .) is very small and can be neglected. In Fig. (l5]-a) we 

■^Conversely, by assuming standard propagation of neutrinos between the source and detector, measure- 
ments of flavor ratios at Earth can be used to set bounds on the initial flavor ratios at the source [26l [47] . 
^a slightly different parametrization We ■ : Wr = 1 '. n : {) has been used in |i47] . 



16 



8.0 
7.5 
7.0 
6.5 
6.0 
5.5 
5.0 
4.5 
4.0 
3.5 
3.0 
2.5 
2.0 
1.5 



— I 1 1 1 1 1 1 r 

No Sterile 

Average on 1 st pair 




1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 
a 



8.0 
7.5 
7.0 
6.5 
6.0 
5.5 
5.0 
4.5 
4.0 
3.5 
3.0 
2.5 
2.0 
1.5 



— I 1 1 1 1 1 r 

No Sterile 

Average on 3rd pair 




1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 
a 

(c) 



8.0 
7.5 
7.0 
6.5 
6.0 
5.5 
5.0 
4.5 
4.0 
3.5 
3.0 
2.5 
2.0 
1.5 



— I 1 r 

Average on 1 st and 3rd pairs 



1 1 1 1 1 r 

No Sterile 




1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 2.3 2.4 2.5 

a 

(e) 



8.0 
7.5 
7.0 
6.5 
6.0 
5.5 
5.0 
4.5 
4.0 
3.5 
3.0 
2.5 
2.0 
1.5 



8.0 
7.5 
7.0 
6.5 
6.0 
5.5 
5.0 
4.5 
4.0 
3.5 
3.0 
2.5 
2.0 
1.5 



8.0 
7.5 
7.0 
6.5 
6.0 
5.5 
5.0 
4.5 
4.0 
3.5 
3.0 
2.5 
2.0 
1.5 



"T" 



"T" 



'no Sterile 
Average on 2nd pair 




2.1 2.2 2.3 2.4 2.5 



— I 1 1 1 1 1 1 r 

No Sterile - 
Average on 1st and 2nd pairs -- 




2.1 2.2 2.3 2.4 2.5 



1 1 1 1 1 1 1 
No Ster 

_ Average on 2nd and 3rd pa 


1 1 

le - 

rs _ 






1 1 
_ 1 1 _ 





1.5 1.6 1.7 1.8 1.9 



2 2.1 2.2 2.3 2.4 2.5 

a 

(f) 



Figure 3: The dependence of R, ratio of ^u-tracks to shower-like events, on a for the stopped-muon 
source with the initial flavor ratio 0:1:0 and power-law spectrum with the spectral index a = 2. 
Assuming the best-fit values for the mixing angles and 6 = 0, the value of R is Rfj_ = 3.07. In 
each figure the red curves represent the case with no sterile neutrino and the green dashed-curves 
correspond to the cases with average on pairs mentioned in the legends. The hatched areas show 
the values that R can take when the input parameters vary in the current uncertainty intervals in 
Table [TJ The two vertical and horizontal dashed-lines show the 10 % and 7 % precisions in the 
measurements of a and R, respectively. 
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Figure 4: The same as Fig. ([3]) except that the input parameters have been varied in the future 
uncertainty intervals of Tabled) Particularly, we assumed 5 % uncertainty interval for the 13-mixing 
angle with the central value near the present upper bound: sin^ ^13 = 0.03. 
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Figure 5: The hatched areas in (a) show the dependence of i? on a for the case of no sterile 
neutrinos (red) and the case with average on the 2nd and 3rd pairs of almost degenerate neutrinos 
(green). The initial flavor ratio Wg '■ : has been varied in the region depicted in (b). In drawing 
this figure we have assumed the best-fit values for the mixing parameters and varied Ag G (0.8, 1). 



have varied n in Eq. f|T5l) such that the separated regions in Fig. ([2]-f) begin to overlap. 
In drawing Fig. ([5]-a) we have varied n G (0.25,0.75). The initial flavor ratios We '■ : 
corresponding to this interval are shown in Fig. ([S]-b). As can be seen, for nearly large 
deviations of the initial flavor ratio from 1 : 1.85 : 0, the hatched areas remain separated, 
which means that the lack of knowledge about the exact value of the initial flavor ratio of 
neutrinos at the source do not affect substantially the capability of neutrino telescopes in 
probing pseudo-Dirac neutrino scenario. 



5 Conclusion 



The new generation of km'^ neutrino telescopes opens a new window to study cosmos via the 
detection of high energy neutrinos predicted to be emitted from the astrophysical objects. 
Because of the extremely large distance of the sources (> Mpc) the flavor oscillation is 
sensitive to the very tiny mass squared differences 10~^^ eV^ < Am^ < 10^^^ eV^. The 
existence of sterile neutrinos with masses almost degenerate with the active ones such that 
the mass squared differences between the active and sterile neutrinos lie in the above region 
is hypothesized in many models (the so-called pseudo-Dirac scenario.) We have studied the 
effect of these sterile neutrinos on the flux of cosmic neutrinos and discussed the capability 
of IceCube to verify the existence of them. In the analysis we have considered different cases 



19 



corresponding to the existence of sterile neutrinos degenerate in mass with different vn such 
that the oscillatory terms driven by these mass squared differences can be averaged out. 

The detection power of IceCube (and other proposed km^ neutrino telescopes), in the 
range of neutrino energies 100 GeV < < 100 TeV, is hmited to distinguishing two types 
of events: /i-track and shower-like events. We have considered the ratio of these events, R 
in Eq. (fT2l) . as the realistic quantity that can be measured in the IceCube. We have studied 
the possibility of using the measured value of i? as a discriminator between the pseudo-Dirac 
scenario and the scenario with no sterile neutrinos. We have considered various sources of 
uncertainties in our analysis. One part of these uncertainties comes from the imprecision 
of the neutrino oscillation experiments in the determination of mixing parameters. For this 
part we considered two sets of uncertainty intervals for the mixing parameters: i) the current 
uncertainty intervals from the performed experiments; and ii) the future uncertainty intervals 
that will be achieved in the forthcoming experiments. Both of these sets have been depicted 
in Table [1] Among the mixing parameters, the measurable quantity R is most sensitive 
to the exact value of the mixing angle ^23- For the uncertainty of this parameter we have 
used the 6 % precision on sin^ 6*23 depicted in Table. [H The other part of uncertainties 
comes from the not completely known mechanism of neutrino production in the sources. 
Many models of neutrino production predict a power-law spectrum for the cosmic neutrinos. 
However, the value of the spectral index in the power-law spectrum depends on the details of 
the neutrino production mechanism and can take values in the interval (1, 3). Also, the ratio 
of the number of electron anti-neutrinos to the number of electron neutrinos is not known. 
We have taken into account all these uncertainties in our analysis. 

The analysis has been done for two different initial flavor composition We '■ : Wr at the 
source: i) pion source 1 : 1.85 : 0; and ii) stopped-muon source 0:1:0. It has been shown 
that with a higher precision of the mixing angles achievable in the forthcoming oscillation 
experiments, in certain cases it is possible to demonstrate or rule out the existence of sterile 
neutrinos hypothesized in the pseudo-Dirac scenario. In these cases the regions corresponding 
to existence and absence of sterile neutrinos are well-separated in the parameter space for 
neutrinos coming from pion sources. For example, it is very promising to probe the case of 
two sterile neutrinos with masses almost degenerate with 1/21 and u^l (see Fig. (El-f))- For 
neutrinos coming from stopped-muon sources, these regions mostly overlap such that their 
discrimination cannot be done without ambiguity. Also, the robustness of these results has 
been tested against the uncertainties in the initial flavor ratio of neutrinos at the source. 
It has been shown that for reasonably large variations of the initial flavor ratio around 
the expected value for pion source, 1 : 1.85 : (see Fig. ([5]) for clarification), the regions 
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corresponding to the existence and the absence of sterile neutrinos remain separated. 
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